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Primordial magnetic fields (PMF) can create polarization B-modes in the cosmic microwave back- 
ground (CMB) through Faraday rotation (FR), leading to non-trivial 2-point and 4-point correlators 
of the CMB temperature and polarization. We discuss the detectability of primordial magnetic fields 
using different correlators and evaluate their relative merits. We have fully accounted for the con- 
tamination by weak lensing, which contributes to the variance, but whose contribution to the 4-point 
correlations is orthogonal to that of FR. We show that a Planck-like experiment can detect scale- 
invariant PMF of nG strength using the FR diagnostic at 90GHz, while realistic future experiments 
at the same frequency can detect 1CP 11 G. Utilizing multiple frequencies will improve on these 
prospects, making FR of CMB a powerful probe of scale-invariant PMF. 

I. INTRODUCTION 

Magnetic fields are prevalent in the cosmic structures around us, in galaxies B ~ 1 /iG with coherence length 
A ~ 1 kpc, in galaxy clusters B ~ 1 — 10 /iG with coherence length A ~ 10 — 100 kpc, and in objects at high redshifts 
z ~ 2 with magnetic field B ~ 10/iG [l[ . Recently there have also been claims of a lower bound on the inter-galactic 
mag netic field, B > 10" 15 G and perhaps a measurement ~ 10 15 G Q, based on the absence of GeV y-ray 

emission in the cascade initiated by TeV y-rays, although it has been argued that plasma instabilities could also 
explain such a non-observation Q. It is possible that these magnetic fields have a common origin from a "seed" 
magnetic field imprinted in the early universe (see [|| for a review) . Magnetic fields may be generated at cosmic phase 
transitions [fjl-fl6j and through specially engineered inflationary mechanisms [It], [HI • 

Detection of primordial magnetic fields (PMF) can lead to important insights into fundamental physics and the early 
universe. Currently, there are upper limits on the strength of PMF from big-bang nucleosynthesis (BBN) [l9l-[23| and 
the Cosmic Microwave Background (CMB) temperature anisotropies [24I27I ]. including their non-Gaussian statistics, 
such as bispectrum and trispectrum [H, [29|. Metric fluctuations induced by PMF are intrinsically non-Gaussian 
because the stress-energy is quadratic in the magnetic field strength B and thus non-Gaussian distributed even if B 
itself is Gaussian. In this paper we study the detectability of PMF through a different observational window, namely, 
the Faraday Rotation (FR) signal they induce in the CMB polarization [3(| ■ The polarization of the CMB field can be 
studied in terms of the parity even E and parity odd B- modes [311433T ] . FR converts some of the primordial i?-modes 
into J3-modes, thus providing a contribution to the -B-mode power spectrum along with the weak gravitational lensing 
and primordial sources, like the actively sought inflationary gravity waves [U, [32|, [34| or cosmic strings (HJ. 

Importantly, spatially dependent FR also couples off-diagonal CMB modes, effectively producing additional non- 
Gaussian signatures in the CMB polarization. The FR induced parity odd correlations (TB) and (EB) must vanish 
in a statistically isotropic universe 1 , where the ensemble average is over many realizations of the stochastic magnetic 
field. However, a particular realization of the FR distortion field that generates a _B-mode from the primordial _E-mode 
will correlate the respective Legendre coefficients E\ m and Bi' m i. In fact, as shown in [38}, it is possible to reconstruct 
the distortion a(n) at a given point n on the sky from specially constructed linear combinations of products Ei m Bii rn i. 
The additional correlations induced by FR also manifest themselves as connected 4-point functions of the CMB, which, 
in turn, provide a measurement of the distortion spectrum C£ Q [39l l4fj| . 

In this paper we discuss the detectability of primordial magnetic fields through estimators based on 4-point correla- 
tions (EBEB) and (TBTB), as well as the 2-point function (BB). In particular, we determine which of the estimators 
has the highest signal to noise for several types of magnetic field spectra and for a range of experimental sensitivities. 
We demonstrate that FR will be a very promising diagnostic of primordial magnetic fields. In particular, future 



This statement is specific to the FR induced parity-odd two-point correlations, which are quadratic in the magnetic field strength B. 
Metric fluctuations sourced by magnetic stress-energy can produce non-zero (TB) and (EB), which arc quartic in B, if the net helicity 
in the magnetic field is non-zero [3a , |37|| . 
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generation of sub-orbital or space-based CMB polarization experiments will be able to detect scale-invariant magnetic 
fields as weak as 10 _11 G based on the measurement at 90 GHz frequency. Measurements at multiple frequencies can 
further significantly improve on these prospects. 



II. FARADAY ROTATION FROM MAGNETIC FIELDS 



Magnetic fields at CMB decoupling will rotate the polarization vector by an angle 
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where f = n e aTd is the differential optical depth, n e is the line of sight free electron density, <tt is the Thomson 
scattering cross-section, a is the scale factor, Ao is the observed wavelength of the radiation, B is the "comoving" 
magnetic field, and dl is the comoving length element along the photon trajectory. We are using Gaussian natural 
units with Ti — c — 1, and the integration limits are from the initial to the final position of the photon. 

Statistically homogeneous and isotropic primordial seed magnetic fields can be generated in the early universe 
during phase transitions [9Hl"tl| , and are described in terms of a two-point correlation function in Fourier space 



(&i(k)&j( k ')> = (2^) 3 <5 (3) (k + k')[(% - kikj)S(k) + iemhAik)] 
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where S(k) and A(k), the symmetric and anti-symmetric magnetic power spectra, are real functions of k — |k|. The 
function A(k) quantifies the amount of magnetic helicity which plays a crucial role in determining the coherence scale 
and the magnitude of magnetic fields as they evolve from an earlier epoch until decoupling. However, only S(k) 
appears in the two-point correlation function of the FR angle, which determines the CMB observables evaluated in 
this paper. As in [4l|, we introduce the dimensionless "FR power spectrum" defined as 



k 3 S(k) 
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This form of the symmetric magnetic spectrum is based on the numerical simulations of causal magnetic field evolution 
in [42|]. In the above, p 7 is the comoving photon energy density, 0b 7 is the magnetic energy density relative to the 
photon energy density, &diss is a dissipation scale above which magnetic fields dissipate, kj is an intermediate inertial 
scale, and 



(4) 



All variables, unless ex plic itly stated, are in comoving coordinates. The exponents n = 5/2 and n' = 3/2 correspond to 
causal magnetic fields [42j , and n w n' sa are expected for magnetic fields generated in an inflationary scenario 
IT^ |. The dissipation scale, fcdiss is not an independent parameter and should, in principle, be dependent on the 
amplitude and the shape of the magnetic fields spectrum. As in [4l| , we assume that k^iss is determined by damping 
into Alfven waves [43|, H and can be related to B c g as 
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where B c g is defined as the effective homogeneous field strength that would have the same total magnetic energy 
density. It is related to f2s 7 via [!l| 



B cff = 3.25 x 10- 6 v^b^ Gauss 



(6) 
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Subsequently, the dissipation scale can be expressed in terms of 51 s 7 as 
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One should be aware of the very approximate nature of the relations ((5|) and ([7]). They are based on the analysis in 
Ref. [lH where small perturbations on top of a homogeneous magnetic field were treated. To extend this analysis to a 
stochastic magnetic field with little power on long wavelengths, Ref. [3 introduced a smoothing procedure and split 
the spectrum into a "homogeneous" part and a "perturbations" part. It is not clear to us if this procedure is valid 
for an arbitrary spectrum, S(k). Instead, we will use Eq. ((5]) as an approximate expression for the dissipation scale. 
We note that this relation also imposes an upper bound on fc/, since ki cannot be greater than fcdiss. 

In practice, both scales {kj and fcdiss) that appear in the definition of the magnetic spectrum ([3]) are likely to be 
smaller than the resolution scale of a realistic CMB experiment. This means that only the large scale tail of the 
spectrum, i.e. the < k < kj range, will be relevant for calculating the shapes of the CMB correlation functions 
on the observable scales. The existence of the intermediate range kj < k < fcdiss an d the exponent n' will affect 
the inferred constraints on f2s 7 only through an overall rescaling. For this reason, in the calculation of the CMB 
correlations, we will set kj = fcdiss, with fcdi ss given by ([7]), and work with a single power law spectrum 
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The bound we obtain on can then be easily converted to the bound on fis 7 via 
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We calculate the FR power spectrum as [4l| 
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2L+V y " 2L + 1 
where 7i(fc) are transfer functions that are independent of the magnetic field 
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Here 77* is the epoch at which the visibility function is maximum, jh are the spherical Bessel functions, and f can be 
easily obtained numerically using public codes such as CMBFAST 45| or CAMB (46j . 



III. THE MODE COUPLING ESTIMATOR AND THE B-MODE SPECTRUM 



FR will rotate the CMB polarization fields generated at last scattering. This introduces coupling between different 
CMB modes which can, in fact, be used to reconstruct the rotation angle map from the observed CMB polarization 
maps [H, H3, S3, 13 . Let T(n), Q(n) and U(n) be the un-rotated CMB temperature field and the two linear polar- 
ization Stokes parameters at angular position n. The temperature is not affected by FR, except for the depolarization 
effect [49l | which would appear as a next order correction and can be ignored in our analysis. Under a rotation of the 
polarization by an angle a(n), the two Stokes parameters transform like a spin two field. Thus, the observed fields 
are 



(Q(n) ± iU(n)) = (Q(n) ± E/(n)) exp(±2ia(n)) 



(13) 
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The observed Stokes parameters can be further combined to form parity-even (_E-mode) and parity-odd (_B-modc) 
combinations which, in the flat sky approximation [3lT[33j are defined by 



[E±iB]{\) = dh[Q(h)±iU(h)]e 



where ip\ = cos 1 (h ■ 1). The relevant 2 ensemble averages of the un- rotated CMB fields can be encapsulated in 

<x(l)>=0, (x*(l)x>(l')) = (2n) 2 6(l-l')Cr' , 



(14) 



(15) 



where x, x' run over the T, E, or B fields, and Cf x are the un-rotated CMB power spectra. 

We would like to isolate the effect of FR on the various correlators that describe the CMB. In particular, gravitational 
lensing causes CMB distortions that can interfere with the computation of the FR effect 3 . The change due to FR 
is the difference between lensed rotated fields and lensed un-rotated fields T, E, B. The change in the CMB fields 
Sx(l) = x(l) — x{l) due to FR is 



5T{\) = 0, 
5B(l) = 2 

SE(l) = - 

= (fi - ipi' 



(16) 



d 2 V 
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E(l') cos 2<pm - B(V) sin 2<p V i a(L), 
B(l') cos 2(p vi + E(V) sin 2<p V i] a(L), 



(27T) 

Thus, due to FR, a mode of wavevector L mixes the polarization modes of 



where L = 1 — 1', and (pn> = 
wavevectors 1 and 1' = 1 L. 

Let us take the ensemble average over multiple realizations of the unrotated CMB fields, while assuming a fixed a 
field. Then, for the rotated variables x ^ x 1 , one can write 



(X*(1K(1')) C MB = / M /(l,l')«(L). 



(17) 



where Jtb = ^Cf E cos 2ip\y , and Jeb = ^\Cf E — C^ B ] cos 2(p\y . Namely, a single realization of a random rotation 
field a will induce parity-odd correlations of types TB and EB that are linearly proportional to the FR angle. If 
we also average over an ensemble of magnetic fields, the above two point functions will vanish for x ^ x' , since the 
expectation value of the magnetic field is zero. 

The mode-coupling rotation (fTT)) imprinted in the CMB by FR imp lies that one can build estimators for recon- 
structing the FR field from the observed CMB. Following Ref. (381 - 4401 148|. l50j . we can define an unbiased estimator 
a xm i(L) for a(L), where x ^ x' , by taking quadratic combinations of different polarization modes weighted by a factor 
F X x'0-i, h)'- 



a xx ,(L) = N X L 
where L = I2 — li , and the normalization 

d 2 h 



d 2 h 



x(\\)x'{L — \{)F xx i{\x,'L - \\) . 



(27T) : 



■f xx ,(h,L-h)F xx ,(h,L-h] 



(18) 



(19) 



is chosen to make the estimator unbiased, i.e. (&(L))cmb = ct(L). The fields x(l) can be obtained from the map of an 
experiment, while the CMB power spectrum of un-rotated but lensed fields can be computed from publicly available 
Boltzmann codes like CMBfast [45| and CAMB (46[. The weights F xx / are determined by minimizing the variance 
subject to the normalization constraint. For xx' = TB and EB the minimization yields 



f x Ah,h) 

h h 



(20) 



2 We assume that the unrotated CMB temperature and polarization fluctuations are Gaussian distributed. 

3 In this paper we do not utilize the frequency dependence of the FR as a tool to differentiate it from other effects. 
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where C xx and Cf x are the observed power spectra including the effects of both the signal and the instrument 

C xx = C xx + A^e' 2 "™™/^ 1 " 2 ) (21) 
where A x is the detector noise and #fwhm is the full- width half-maximum (FWHM) resolution of the Gaussian beam 

The variance of the estimator can be calculated as 
Var(<w(L)) = (a xx ,(L)a* xx ,(V)) 

= NfNff' J -^L J ^(x^x'iL-^x^x'iL' -l 2 ))F xx ,(h,L-h)F xx ,(h,L' 

= (2Tr) 2 5(L-V){CZ a +NZ x '}. (22) 

In the last line, the first term is the desired FR power spectrum and the second term is the noise - also referred to as 
the Gaussian noise - in the reconstruction of FR. Note that the noise turns out to coincide with the normalization of 
the minimum variance estimator as given by Eq. (fl~9| [38l . IZol ItR [48|, and is independent of FR. In principle, there 
are higher order noise terms (referred to as non-Gaussian noise terms) which depend on FR, however, these terms 
are sub-dominant in comparison to the Gaussian noise. The signal-to-noisc for detecting a spatially varying FR angle 
a(n) using the estimator (|T8|) is given by [H, |39| 

^(2* + l)N£ir) . ( 23 ) 



N J , ^ 2 v J \N, 



1=2 



where Cf a is the rotation angle power spectrum, and we ignore the contribution of Cf a to the variance (|22|) as it is 
negligible compared to N xx . 

We also consider the case when x = x' = B, and the the resultant BB correlation is quadratic in the rotation field. 
Averaging over the ensemble of magnetic fields gives the FR induced CMB B-mode power spectrum: 



n BB A f d 1' rma/i££ „„ c , 2 TO,„ 1 



2V 

(24) 
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The signal-to-noise in this case (accounting for the i?-modes from weak lensing but assuming no contribution from 
inflationary gravity waves) is given by 

where Af t BB ~ Cf ' ensms + (jB B nolsc an( j |- ne no j se power spectrum is (j BB ' nolsc — Ap exp(/ 2 8p WHM /81n2), where 
Ap is the instrument noise for polarization. 

One might wonder if gravitational lensing, which also generates off-diagonal correlations in the CMB with a leading 
order contribution to the trispectrum, might bias the mode coupling FR estimator (|18[) . However, it has been shown 
that the effects of gravitational lensing and rotation are orthogonal and hence lensing does not bias estimates of 
rotation (39l |40| . Lensing does, however, increase the variance of the estimator, i.e. the Gaussian noise given by 
Eq. (fT9|) is enhanced. For example, the noise of the EB estimator, N EB , depends on the observed CMB E and B 



mode power spectra, and is an integration over {Cf E +N BE )(C BB +N BB ). Further, the f?-modes due to lensing scale 
as white noise up to I ~ 1000, and correspond to a noise level of 5/iK-arcmin. Therefore, the increase in the variance 
is sub-dominant for experiments with Ap > 5/iK-arcmin. However, for small instrumental noise, Ap < 5/LtK-arcmin, 
lensing B-modes become important, saturating the variance to ~ 10~ 6 deg 2 even for an ideal experiment. 

In the next section, we compare detectability of the FR signal in upcoming and future CMB maps with the mode- 
coupling estimators (|18|) vs. using the B-mode spectrum (|24[) . Interestingly, we find that one method can out-perform 
the other depending on the value of the magnetic spectral index. 



IV. DETECTABILITY OF PRIMORDIAL MAGNETIC FIELDS 



To forecast the detectability of FR we consider three experimental setups: a Planck-like satellite [52[ (El), a ground- 
or balloon-based experiment realistically achievable in the next decade (E2), and a future dedicated CMB polarization 
satellite (E3). We will adopt the following parameters for the three experiments: 
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FIG. 1: The noise L(L + 1)N£ X /2tt as given by Eq. (fT9l) for the EB (upper panel) and TB (lower panel) estimators as a 
function of multipole L for three experimental setups (El, E2, E3). Also plotted is the FR power spectrum L(L + l)C^ a /2n 
for the causal magnetic spectrum with 2n = 5 and for a nearly scale invariant spectrum with 2n — 0.1. The amplitude of the 
FR spectrum is set by the observed frequency of vo = 90 GHz and Qbj — 10~ 3 for the causal case (2n = 5), and &b\ — 10 -4 
for the nearly scale invariant case (2n = 0.1). 



El: a Planck-like satellite |52j with noise level Ap = 60/iA'-arcmin and 6fwhm = 7' 
• E2: a next decade sub-orbital experiment with Ap = 3.0/xif-arcmin and 6fwhm = 1' 



E3: a CMBPol-like instrument with A P = ^/2^K -arcmin and 6fwhm = 4', typical of the proposed future 
space-based CMB experiments [H, H3] 



The forecasts directly depend on the fraction of the sky covered by the experiment, / s ky, which is close to unity for 
El and E3, and will be smaller for E2. We quote our bounds subject to specifying / s k y which only appears under a 
quartic root in the bounds on B c s . 

In Fig. [1] we plot L[L + 1)A^|; X /2w vs multipole L with the noise for the EB and TB estimators given by Eq. p^|) 
for experiments El, E2 and E3, along with the FR spectrum L(L + l)C£ Q /27r for two choices of the magnetic spectral 
index 2n. The amplitudes of the rotation spectra are normalized to f2p 7 = 10~ 3 for the causal case (2n = 5) and 
^B-y = 1CP 4 for the nearly scale invariant case (2n = 0.1), with f2p 7 related to J7b 7 via Eq. (QJjJ). For all the estimates 
in this paper we adopted an observational frequency of vq = 90 GHz, but one can easily scale the signal to other 
frequencies using Eq. ©. For all three experiments under consideration, the EB estimator is more sensitive than the 
TB, with the noise Nff° staying roughly constant up to L ~ 1000. Although we do not show L — in the plot, these 
estimators can also be used to estimate the detectability of uniform rotation. 

In Fig. [2] we plot the noise contribution to the variance of the -B-mode spectrum, along with the FR induced 
B-mode spectrum for J7p 7 = 10~ 3 for 2n = 5, and Vls-y = 10~ 4 for 2n = 0.1. In the latter case, the shape of 
the i3-mode spectrum is a close copy of the underlying E-mode, except that the FR induced spectrum falls off as 



L(L + l)Cff ex L 1 at high L |4l|,l55| compared to the exponential fall off of the primordial _E-mode. The L 



2)i- 



tail implies a sharply rising spectrum for the causal case with (2n = 5), with most of the power concentrated near the 
dissipation scale Ldiss ~ 10 4 Mpc x kdh 
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FIG. 2: The B-mode polarization pixel noise, L(L+1)C"" ' nolsc /2n, for experiments El, E2, and E3 as a function of multipole 
L. Also shown is the -B-mode spectrum induced by FR for a causal magnetic spectrum with 2n = 5 and for a nearly scale 
invariant spectrum with 2n = 0.1. The observed frequency is set to vq — 90 GHz, Qb-, = 10~ 3 for the causal case, and 
Q,Bj = 10~ 4 for the scale invariant case. 



To forecast the minimum detectable magnetic field energy fraction, we define it as f2p 7 for which S/N in Eqs. (|23[) 
and (|25p is unity. Note that 0p 7 determines the dissipation scale fcdi SS via Eq. (|7J|. We perform this forecast for 
each estimator, for different experiments, and for several choices of 2n. We restrict the maximum multipole to 
L max = 10000. Fig. [3] shows the minimum detectable fip 7 as a function of the instrument noise Ap for a fixed beam 
©fwhm = 8' for two choices of 2n. As one can see, at high noise levels (Ap > 100/iK-arcmin), the B-mode power 
spectrum tends to cither give comparable or better constraints than the EB estimator. In particular, it is always the 
better probe of the causal primordial magnetic fields. However, for upcoming polarization sensitive experiments with 
lower levels of noise, the EB estimator will become almost comparable to the B-mode power spectrum for causal 
fields, and outperform it when probing scale-invariant FR fields. 

In Fig. |]we plot the minimum 51p 7 that will be detectable by the El, E2 and E3 experiments depending on the 
value of the magnetic spectral index 2n. We see that for all experiments, the EB estimator begins to outperform the 
B-mode spectrum when 2?i < 1, while the TB estimator is always the third best. 

The relative strengths of the three estimators, demonstrated in Figs. [3] and 21 can be understood as follows. 
Generally, the EB and TB estimators have a larger number of independent modes contributing to the signal than the 
-B-mode spectrum. Thus, in principle, it is not surprising if they result in a higher signal to noise. However, whether 
that is the case depends on the experimental noise level, and the distribution of power in the given combination of 
CMB fields and in the magnetic field. For a scale-invariant PMF spectrum, the -B-mode is essentially a copy of the 
-E-mode, with most of the B-mode power being on scales where the -E-modes are also most prominent. This results 
in a strong correlation between E and B for scale-invariant fields. In the case of the TB correlation, the underlying 
T and E (B is obtained by a scale-invariant rotation of E) fields peak on rather different scales. Namely, T peaks at 
£ ~ 200 while E peaks at i ~ 1000. In other words, the intrinsic correlation between T and E is already suboptimal, 
translating into a lesser correlation between T and B. Thus, for experiments with a sufficiently low noise Ap, such 
as El, E2 and E3 considered in this paper, the EB estimator performs better than TB for scale- invariant fields. This 
would not necessarily remain true if polarization measurements had a significantly higher experimental noise. 
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FIG. 3: Comparison of the three estimators, BB (solid red), EB (dashed green) and TB (dotted blue). Plotted is the minimum 
detectable magnetic field amplitude f2.e 7 as a function of the experimental noise, with 9fwhm = 8'. The upper panel is for 
the nearly scale invariant case with 2n — 0.1, while the lower panel is for the causal case with 2n — 5. 



For the blue causal spectra, the FR power is concentrated on very small scales, far away from the scales at which 
any of the unrotated CMB fields have significant power. This means that the -B-modes in the observable range are 
obtained either by the rotation of -©-modes far away from their peak power scale, or by a rotation of peak iJ-modc 
by a negligible angle. This means that E and B fields peak at very different scales, with their correlation being close 
to zero over the observable scales. In this case, we see that the S-mode spectrum, i.e. the BB correlation, has the 
highest signal to noise. 

When interpreting the forecasted bounds on the magnetic field energy fraction or the effective magnetic field strength 
in Figs. [3] and IH several points must be kept in mind: 

1. The constraints are on Qb*/, obtained after setting fc/ = fcdiss, with the dissipation scale determined from Eq. (|7|). 
For scale-invariant fields there is no difference between f2s 7 and ils-y, since the factor relating them in Eq. (|T0|) 
goes to unity when 2n — > 0. Also, for scale-invariant fields, the effective field B c g defined via Eq. ^ is the same 
as the commonly used B\, which is the field smoothed on a given scale A. Thus, our forecasts of the minimum 
detectable B c r for scale-invariant fields can be directly compared to most other bounds in the literature. 

2. For causal fields, the bound on flsj will generally overestimate the magnetic energy fraction, since it assumes 
that the spectrum will keep rising at the same steep rate (2n = 5) all the way to the dissipation scale, which is 
much smaller than the smallest scale directly probed by CMB experiments. Simulations [42[ suggest that the 
spectrum must become less steep, with 2n' = 3 in the range fc/ < fc < fcdiss, implying a smaller net magnetic 
energy fraction ils-y- Since the value of fc/ is not well-known at this point, we chose to quote our bounds in 
terms of ri_g 7 , while keeping in mind that bounds on f2s 7 for causal fields will generally be tighter. 

3. Strictly speaking, our bounds are on the fraction in magnetic fields at the time when the initial conditions for the 
transfer functions (|12p were set, which is a time close to last scattering. While it is expected that the magnetic 
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FIG. 4: The minimum detectable magnetic field amplitude f2s 7 as a function of the magnetic spectral index 2n for the three 
estimators, BB (solid red), EB (dashed green) and TB (dotted blue). The top panel is for El, the middle panel is for E2 and 
the lower panel is for E3. 

fields arc effectively frozen-in between the BBN and last scattering, with a relatively slow time evolution of the 
dissipation scale, this is still an approximation. 

4. The bounds are based on using a single frequency band. Using several bands will improve the constraints. 

Generally, CMB is not very sensitive to magnetic fields with blue spectra because most of the anisotropics are 
concentrated on very small scales. This is what Figs. [3] and 2] are showing too. However, looking for the FR signatures 
at many frequencies can potentially improve the existing CMB bounds on causal fields by a large factor. We leave 
this question as a topic for future exploration. 

In the case of scale-invariant fields, current bounds on the magnetic field strength from WMAP are at a level of a 
few nG 0, [H, I27l - l29| . These bounds are based on the anisotropics induced by the metric fluctuations sourced by 
magnetic fields, and ignore the FR effect. In Refs. [26|, |4l[ the WMAP bound using FR was obtained at the 10~ 7 G 
level. As one can see from Fig. [4j Planck (El) can almost match today's bounds for scale invariant (n = 0) fields 
using the EB estimator at only one frequency, while future probes, such as E2 and E3, can improve the bounds by 
two orders of magnitude! This suggests that the mode coupling estimators of FR will be a very powerful, if not the 
most powerful, direct probe of scale-invariant magnetic fields at the time of last scattering. 

V. SUMMARY 

A primordial magnetic field present at and just after last scattering will Faraday-rotate the plane of polarization of 
the CMB photons. The FR will create i?-modc polarization even in the absence of primordial sources, such as gravity 
waves. In addition to a B-modc autocorrelation, FR also couples different modes of the CMB fields, generating specific 
off-diagonal correlators. Estimators of the polarization rotation angle, such as (|18|) . can utilize these non-Gaussian 
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features to reconstruct the FR map. One can construct four such estimators containing products of two CMB fields, 
one of which contains polarization: TE, EE,TB, and EB. Of these four, the first two receive a large contribution to 
their variance from the usual scalar adiabatic Gaussian perturbations which makes it harder to find the FR signal. 
In this paper, we have considered the last two and found that the EB estimator has the highest signal-to-noise. 

For causal magnetic fields, which tend to have very blue power spectra, the -B-mode power spectrum has a higher 
signal-to-noise due to reasons explained in the previous section. However, the EB estimator performs better for 
scale-invariant fields. In addition, there are certain advantages in using mode-coupling based estimators, such as EB, 
over the traditional .B-modes power spectrum. For instance, there are other sources that can generate £?-modcs, 
such as weak lensing, patchy reionization, inflationary tensor perturbations or cosmic strings, and, in fact, metric 
perturbations induced by the magnetic fields. Hence, one has to separate the FR induced B-modes either by using 
the frequency dependence of FR or features in the B-mode spectrum. Interestingly, while patchy reionization and 
lensing also generate off-diagonal correlations in CMB, their contributions are "orthogonal" to the features imprinted 
by FR |48| . Hence, mode-coupling estimators do not suffer from contamination from other contributions. This means 
that a larger part of the information in the frequency dependence can be used for systematic checks and to separate 
from other foreground contamination. In addition, mode-coupling estimators can be used to reconstruct the map 
of FR. This FR map can be used to cross-correlate with other tracers of magnetic fields, including CMB maps and 
surveys of large scale structure. Such correlations studies are useful for systematic checks and for increasing the 
signal-to-noise. 

We have found that a Planck-like experiment at 90GHz can detect scale-invariant PMF of a few nG strength, which 
is comparable to the ~nG sensitivity forecasted for Planck based on information in the CMB temperature anisotropies 
p7j . Future CMB experiments will be able detect scale-invariant fields as weak as 10~ n G at 90GHz. Thus, FR 
should become a leading diagnostic of PMF when analyzing future CMB polarization data. 
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